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ABSTRACT 


This thesis extends the development of the Generalized Value System (GVS), used 
in the Airland Advanced Research Model (ALARM), as an on-going research effort at 
the Naval Postgraduate School. Specifically, the problem of determining the multidi- 
mensional mapping of the state variables that represent the condition of an entity into 
the power function is addressed. The methods described in this thesis provide a means 
of aquiring this mapping function by the use of a Degraded Power Function (DPF). 
The DPF provides a basis for estimating the future state of an entity based on the state 
(condition) of the entity, virtually eliminating the exponential decay function presently 


In use. 


THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While every effort has been made, within the 
time available, to ensure that the programs are free of computational and logic errors, 
they cannot be considered validated. Anv application of these programs without addi- 
tional verification is at the risk of the user. 
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I INTRODUCTION 


Modeling of the decision-making process in present combat simulations continues 
to be the driving factor behind the improvements in such models as the AirLand Ad- 
vanced Research Model (ALARM), being developed at the Naval Postgraduate School. 
ALARM is designed to evaluate a U.S. Corps versus a Soviet Front. The decision- 
making process unquestionably impacts upon the accuracy of any combat simulation. 
It is therefore, very important to insure the decision-making process takes into account 
all factors impacting on the ultimate goal: mission accomplishment. In order to have 
realistic simulations, there must exist a reliable data base with all the information needed 
to run the simulation. Obvious shortcomings exist if the data base 1s not representative 
of the “real world”. 

To determine the combat effectiveness of a unit and make any conclusions or rec- 
ommendations regarding its ability to accomplish the assigned mission requires the 
decision-maker to answer several questions: 

@ What is the assigned mission ? 

© What resources are available ? 

e What is the enemy situation ? 

® How much time ts available ? 

e What is the condition of personnel ? 

e What type of terrain will be encountered ? 

e What are the capabilities of the friendly units ? 


e What are the enemy capabilities ? 


The focal point of the decision-making process is the decision-inaker’s appraisal of his 
units’ combat effectiveness (ability to accomplish its assigned mission), based on the 
answers to the questions given above. In general, the problem is to describe the relative 
importance of the various resources available and decide if the unit can accomplish its 
assigned mission. Staff planning must take into account how the available resources 
constrain the combat operation. Once the decision maker has gathered both the infor- 


mation needed and the recommendations of his staff, he then determunes the best course 


of action to take in order to accomplish the assigned mission. The more information the 
decision-maker has at his disposal the better he can make a critical decision with rea- 
sonable certainty. 

This thesis presents a dynamic approach of mapping the available resources (state 
variables) into the unit’s power. In order to accomplish this mapping function and in- 
corporate the state variables in the computation of power there ts first a need to estab- 
lish a data base which relates the state variables to combat effectiveness. 

There are several means of determining a unit’s effectiveness. Undoubtedly, the 
best method is actual combat. However, there are other methods, short of warfare, 
Which provide a viable determination of combat effectiveness. One of the better metli- 
ods is based on the judgement of expert deeision-inakers. 

The primary goal of this thesis is the development of a usable and realistic model 
which will describe the relationship between the state variables, representing the re- 
sources available to a combat unit, and the resulting combat effectiveness of that unit 
and the unit’s abilitv to accomplish its assigned coibat mission. It is important to note 
that the basic methodology outlined tn this paper 1s apphcable to a wide variety of situ- 
ations in which the decision-maker 1s required to view the current situation and make 
an assessment of the unit’s combat effectiveness. 

Background information is presented in Chapter II, describing the techniques and 
systems being used in ALARM. Chapter III discusses different techniques of gathering 
information from judges, along with a detailed explanation of the teclinique used in tls 
thesis. A substantial effort was made to determine the best techniques of asking the 
right questions and the rnethod chosen was deenied the best of the methods investigated. 
In Chapter IV the data is analyzed and the best fitting model is discussed, along with 
an example of the sensitivity of tlie model developed. Finally, future applications are 


suggested, based on the methodology developed in this thesis. 


I. BACKGROUND 


In present day military operations there has emerged a term called Airland battle, 
which is becoming the doctrine of the U.S. Army. This doctrine states that battles will 
not be limited to the forward fighting elements but will involve the total battle zone to 
include the rear areas. In order to more closely examine this doctrine and its effects on 
battle outcomes, computer simulations are used to examine the effects of this doctrine. 
This need for evaluation motivated the design of the AirLand Advanced Research Model 
(ALARM). ALARM, under development at the Naval Postgraduate Scliool since Jan- 
uary 1985, focuses on investigations of the effectiveness of the Airland battle doctrine 
for a U.S. Corps battle. [Ref. 1] ALARM 1s capable of man-in-the-loop decisions but 
is primarily being developed as a systemic (1.e., no man-in-the-loop) U.S. Corps versus 
Soviet Front model. The ALARMI concept incorporates three basic methodologies: 
Time Domain Networks, Cartesian Space Networks, and Generalized Value System. 

The function of time domain networks is “to handle the planning function or ac- 
tivity Within the ALARM system. The time domain network consists of nodes and arcs 
(or links) connecting the nodes. The key to time domain networks ts that arcs do not 
represent distance, but represent the passage of trme or the completion of a sub-activity 
which leads to the completion of the entire activity represented by the network.” [Ref. 
1: pp.6-8] 

Cartesian space networks consist of “a series of networks each representing a dif- 
ferent aspect of the battlefield. Each network represents physical connections between 
points on the battlefield. ALARM will have three or more Cartesian space networks.” 
[Ref. I: p.8] The three networks that have been tdentified are as follows: 

1. Terrain and Transportation networks [Ref. 2] 


2. Communications network 


3. Logistics Resupply network 


The Generalized Value System (GVS), an axiomatic value system used as a tool for 
evaluating the power of entities at present and future times, was first proposed m 1985 


by Professor Schoenstadt of the Naval Postgraduate School [Ref. 3] and was expanded 


in 1986 by Robert Kilmer [Ref. 4]. This procedure provides a means of forecasting fu- 
ture states of entities in continuous time, allows different entities to be compared, and 
works for all entities on the battlefield. Essentially. “GVS is a procedure for quantifving 
the capabilities and importance of entities on the battlefield at some future time (t +x). 
It does this by using algorithms to predict future entity or situation states at time (t+ x) 
based on the situation at time (t). Thus GVS provides the framework for forcasting fu- 
ture states of entities in continuous time.” [Ref. 1: p.8] 

Inherent in this procedure is the concept of umit effectiveness and assigning that 
unit (entity) a value which 1s consistent with other units (entrtres) on the battlefield. This 
value, called “Power,” is the value that GVS deterinines based on the usefulness of an 


entity at that time. 


A. POWER 
Power is a subjective value given to an element or unit on the battlefield. In GVS 
it combines both inherent and derrved power. This value is then used as a planning 
factor in the combat simulation. Accuracy m assigning values to elements is crucial to 
the reliability of this methodology. IIence, a great deal of time and effort has been spent 
in the assignment of such values. Several of the techniques used in determining values 
are firepower scores, user input values, and Generalized Value Svstem (GVS). Firepower 
scores use an aggregated approach to determine the value of an entity (Le., the value 
of an entity depends on the additive value of all of its weapon svstems). User Input 
Values are assigned along with priorities of target selection bv the user. Generalized 
Value Svstem (GVS) assigns a value to different entities on the battlefield, to inelude 
non-combat entities such as combat support units, bridges, etc... Each value 1s compa- 
rable to other values. In GVS an entity relates to a unit or even a bridge and can have 
two types of power: inherent and/or derived. A combat unrt will normally have a com- 
bination of inherent and derived power, while a combat support unit may only have 

derived power. 

|. Inherent power 

Inherent power is the ability of a unit (entity) to affect the outcome ofa battle 
by the disruption, destruction, and confusion of tlle enemy forces. I[t also includes the 
effectiveness and survivability of a unit. Kilmer defines inherent power as “the ability 
to directly affect the states of enemy entities or of entities that the enemy 1s using or 
planning to use (e.g., a bridge).” (Ref. 4: p.32 ] Inherent power tvpes that need to be 


defined in order to focus the direction of this paper are Basic Inherent Power (BIP) and 


Situational Inherent Power (SIP). Basic inherent power (BIP) is the inherent power that 
a unit possesses when it is at full strength, in position and ready to perform its most 
likely mission. Situational Inherent Power (SIP) 1s the inherent power that a unit 1s 
predicted to have at some future time, t. 

To calculate the Situational Inherent Power (SIP) of an entity, the time the 
entily is expected to be ready to accomplish its mission is estimated from the time do- 
main networks. This time is called the available tune (7,). Other aspects that must be 
determined are the rate of attaining readiness (D), the present time (c,), and the Predicted 
Adjusted Basic Inherent Power (PABIP). Once these values are known a unit's SIP can 


be computed using equation 2.1: 


SIP(t|t,) = PABIP(SX1(t,)) xe "4-O for OS tS ty, (2.1) 
Ser |t,) = PABIP(¢|t,) onal ts 
moliere: 


D is the rate at which an entity is attaining readiness, 
t, = Time the entity will be available to execute the assigned nussion, 


meee Present tine, and 


te 


SN1(c) = State (condition) of entity X1 at time t expressed as a vector of state valucs. 


It can be seen that the SIP of an entity is dependent on the PABIP, but the PABIP is 
dependent on the Adjusted Basic Inherent Power (ABIP). The ABIP is the percent de- 
gradation at any time, t. Therefore ABIP = (percent degradataion) x BIP. From this 


the PABIP can be calculated using equation 2.2: 


PABIP(XA( t )|SXM(t,)) = ABIP(SX\ | eee es on <6 =. 


rN 
rN 
ty 

wane? 


where: 


L is a rate of change which may have several components such as attrition, re- 
enforcements, and logistics effects, and 


f, 1§ the end of the planning horizon. 


This equation uses an exponential decay function to calculate PABIP of an entity over 


the entire planning horizon. Note that, in general, ¢, may occur either before or after 


te. This thesis proposes a inethod of determining the PABIP based on the state of an 
entity and not on an exponential decay function. This method will be discussed in detail 
in Chapter [II. 
2. Derived Power 

Derived power consists of the power a unit has when it supports another unit 
or provides the services needed by a unit in accomphshing its assigned mission. Kilmer 
defines derived power of an entity as “the power it possesses because of its ability to in- 
fluence the states of other friendly entities or of entities that its forces are planning to 
USC. qtincl. 4: poo 

The estimation of power in modern computer wargaming exercises has been 
evaluated and questioned for many years. There is a great need to have the best rep- 
resentation of a units’ power. In the present ALARM the power is based on the value 
GVS gives the entity and GVS uses an exponential decay function which 1s not directly 
related to the present state of an entity. ALARM does have formulations to determine 
the units’ present strengths and thus the the state of the unit (entity). The state of the 
entity takes into account the factors impacting on the ability of the entitv to accomplish 


its assigned mussion. 


B. THESIS MOTIVATION 

In order to inake combat simulations more accurate and have less arbitrary as- 
signment of values for power, there 1s a need for a multivariate mapping function of state 
vectors to unit power. This mapping function will allow a future state decision process 
to replace the exponential function presently being used in deterinining the PABIP. 
Conversely, this mapping function will allow GVS to determine the future state of an 
entity based on the state variables which are representative of the entity. This mapping 
function is discussed in Chapter III and is based on the collective judgement of a large 
group of decision-makers vice an arbitrary assignment of values. The coinbined expert- 
ence of the decision-makers to determine the future state of an entity through the map- 
ping function creates a more accurate combat simulation and more accurate results. 

Power, therefore, should not be an arbitrary determination made by one individual 
or a sinall group of “experts,” but should be decided by a large group of actual 
decision-makers. Since there is no data base established from which to determine the 


power of an entity, there exists a need to establish such a data base. This need to more 


accurately represent the actual battlefield by developing a multivariate mapping function 
of state variables into unit power. This will enable a better determination of the future 
state of an entity and make GVS a more usable technique. Thus is the motivation for 


tis thesis. 


I. THE EXPERIMENT 


As stated in Chapter I, the primary motivation for this thesis is to develop a 
multivariate mapping function of state variables into unit power in order to better de- 
termine the future state of an entity. In conducting an experiment to develop tlus 
function it was first necessary to determine a technique that would allow for the rnost 
precise, reliable responses from a select group of judges. First a discussion of the de- 
termination of combat effectiveness by a decision-maker is needed. 

To determine the combat effectiveness of a unit and draw any conclusions regarding 
Whether or not the unit can accomplish its assigned mission requires the decision-inaker 
to look at several questions: 

e What is the assigned mission ? 

@ What resources are available ? 

@ What is the enemy situation ? 

¢ flow much time is available ? 

e What is the condition of personnel ? 

e What type of terrain is being defended ? 
e What are the capabilities of the unit ? 


e What are the enemy capabilities ? 


As pointed out in Chapter I, the focal point of the decisiton-niaking process is the deci- 
sion makers’ appraisal of his units’ combat effectiveness (ability to accomplish its as- 
signed mission) based on the answers to the questions listed above. The capabilities of 
the friendly units are critical in this planning process (available resources can determine 
the feasibility of any defensive operation). Once the decision maker has gathered both 
the needed information and the recommendations of his staff, he then determines the 
best course of action to take in order to accoinphish the assigned mission. The more 
information the decision-maker has at his disposal, the better he can make a critical de- 
cision with reasonable certainty. 

The questions posed above are a key element in the decision-making process and 
as such they are the descriptors of the unit (1.e., they tell the decision maker the STATE 


of his unit). These descriptors will be called the state variables for this study and the 


STATE of the unit will be based on the value of each variable. In order to better uu- 
derstand this concept a thorough description of the systern and state variables 1s re- 
quired. It is important to provide as detailed and comprehensive a list of variables 
linpacting on coinbat effectiveness as feasible, so that any simplifving assumptious and 
State variable onussions made during the course of problem definition and analvsis inay 
be directly related to the initial svstem description. 

In general the problem is to describe the relative importance of the various state 
variables in a subjective decision process. It is liypothesized that the decision maker 
considers his units’ capabilities and his vulnerabilities in deciding on his units’ ability to 
accomplish its mission. There are 7 variables that are considered, in this study, to aflect 
combat outcome. 

1. Personnel strength (foxhole strength). 
Available ammunition. 
Key operational weapon systems. 
Fuel and oil products available (POL). 
Chain of command (key leaders). 


Level of experience of unit and leaders. 


eo a ee 


Enemy force capabilities. 


It would be ideal to be able to develop this model utilizing all of these variables and 
evaluate them all at many levels, but that would clearly be unmanageable. We will 
therefore reduce the number of these variables to a manageable level. The mission of 
the unit will be limited to deliberate defense. The chain of command will not be evalu- 
ated since when evaluating personnel strength it is assumed that this wil! include kev 
leaders. As experience of the unit and of the leaders is not quantifiable, it will therefore, 
not be evaluated. Enemy forces will be a Motorized rifle regiment, and therefore will 
not change. This leaves four variables which will be used in this model: 


1. Current personnel strength, as a percent of the authorized Table of Organization 
and Equipment (TO&E), (The TO&E is a document that reflects all of the types 
and quantities of personnel and equipmeut a particular unit is authorized), 


2. Current ammunition available, as a perceut of the authorized amount, 
3. Current combat vehicles operational, as a percent of the authorized TOE, 


4. Current POL fuel and oil products available, as a percent of the authorized 
amount. 


It must now be determined at what levels these variables will be evaluated. Based 
on discussions with combat arnis officers and other experts, as well as review of an ear- 
lier research paper bv Elizabeth Etheridge and Michael Anderson [Ref. 5], it was deter- 
inined that personnel and combat vehicles would be evaluated at three levels (100%, 
75%, 50%) and ammunition and POL would be evaluated at four levels (100%, 75%, 
50%, 25%), constituting a factoral design. A factoral experiment is one in which all 
levels of a given variable are combined with all levels of every other variable in the ex- 
periment. A 4* factoral experiment 1s a case where k factors are evaluated each at 4 
levels. It can easily be seen that a 4* factoral experiment requires 4* = 256 separate 
unique combinations, each representing a different possible value of combat effective- 
ness. In the case of the model being developed where two variables are being evaluated 
at three levels and two vaniables are being evaluated at four levels, there exists a 


(3?) x (47) factoral experiment requiring [44 separate evaluations. 


A. STATEMENT OF THE PROBLEM 

Deterinining the relative unportance of the various system state variables in a 
decision-maker’s estimate of particular variables relative to his units’ ability to accom- 
plish its assigned mission is the ultimate objective of this study. In determining the rel- 
ative value of these variables, a subjective evaluation was used in order to establish a 
data base. Several techniques were considcred in the establishrnent of this data base: a 
continuous response scale, the paired comparison test, ordinal judgements and 


categorical judgements [Ref. 6]. Each 1s discussed briefly below. 


B. CONTINUOUS RESPONSE SCALE 

A continuous response scale allows the judges to respond on a continuous scale 
from 0 to 100, based on the judges’ feelings of the items he is evaluating. This type o! 
test could lead to an Analvsis of Variance (ANOVA) design which is commonly used to 
determine interrelationships. The problem with this design is that judges may find it 
difficult to make a judgement on a continuous scale froin 0 to 100 and for this reason 


this method was rejected. 


C. PAIRED COMPARISON TEST 

In the paired comparison test, the judges are presented with pairs of items being 
evaluated (instances) and are then asked to determine which possesses the greater value. 
lor a paired comparison test the judges are asked to compare each of the possible pairs, 


which involves n(n-1)/2 comparisons. For this study it would involve 144(144-1)/2 = 


10 


10,296 comparisons. This test allows each variable to be compared with all other vari- 
ables independently helping to determine any miterrelationships between variables. The 
problem with this method, however, is that it requires a large number of comparisons. 
Since the variables are paired together, a higher order interrelationship cannot be deter- 


mined. Therefore, it was not used tn this study. 


D. ORDINAL JUDGEMENTS 

In making ordinal judgements, judges are asked to rank instances based on their 
“feelmgs” about an instance and subsequently set up a rank for all of the instances that 
they are evaluating. The problem with this technique is that a judge may not be able to 
discriminate between two instances in terms of the value he mav feel they possess, and 
when comparing a large number of instances the problem ts confounded. Also, some 
judges may not rank all mstances. This technique was not used since the number of in- 


stances each judge would have to rank would be prohibitive. 


E. CATEGORICAL JUDGEMENTS 

This method requires the judges to select the category that thev think best repres- 
ents an instance. The categories are assumed to be a mutually exclusive set of successive 
intervals on the variable’s scale. There are descriptors with each category that help the 
judges with their task, (1.e., strongly agree, agree, no opinion, disagree, strongly disa- 
gree). Normality is an assumption of this technique, where the judges feelings about the 
scale value of unit effectiveness for an mstance is normally distributed, as seen in Figure 
L 

The categorical judgement technique did not require an exhaustive amount of work 
by the judges and discussions with several experts in the field at the Naval Postgraduate 
School provided greater assurance that choosing this technique would result in the final 


development of a usable model. 


F. CONSTRUCTING INTERVAL SCALES FROM CATEGORICAL RESPONSES 

Using the categorical judgements technique to transform the judgemental responses 
to an interval scale was determined to be the best technique. Hence finding an mterval 
scale based on the responses of the judges requires ten steps. [Ref. 7: p.6] These steps 
are summorized here and given in detail, in Section G below. This technique also re- 


quires 4 assumptions. 
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Figure 1. Normal Population of Judges Feelings of Unit cffectiveness 


1. The judges’ “feelings” about the scale value of scenario (instance) 1 are assumed 
to be a normally distributed random variable with mean mw, and variance o?. 
2. Categories are a mutually exclusive set of of successive intervals which collectively 


exhaust the continuum of the scale. 


3. The judges feelings about a categorv’s upper bound as a normally distributed 
random variable, so for category } the upper bound would be a normially distrib- 
uted random variable with mean wv; and variance o} . 


te 


All category bounds are assumed to have the same variance, so for allj, a? = c. 


The respouses are recorded in a table that lists all of the raw frequencies for each 
instance. The raw frequencies are then divided by the total number of judges who an- 


swered that instance to determine relative frequencies. The relative frequencies are then 


used to calculate the cumulative relative frequencies which are likened to the cumulative 
frequencies under a normal curve. The cumulative relative frequencies constitute the P 
array, Where P = || p;,, ||,and p,, is the cumulative relative frequency for instance 1, 
and category j. These p,, are then grouped with other instances having similar category 
values (i.e., instances with all values in the categories A, B, and C are grouped with all 
other instances with values in categories A, B, and C). Once the instances are grouped, 
all values of p,,2 0.98 and p,,< 0.02 are renioved to avoid anv undue influence by a 
small number of judges. Treating these values of p,, as leftward areas under a Normal 
Mempecurve, a table of the Normal CDF values ts developed. These values are the Z 
array used in the computations, where Z = |] z,, ||, and z,, is the value from the normal 
table for instance 1, and category j. The row, column and grand average are calculated 
from this Z array, along with the estimate of the standard deviation. The scale values 
of the instances are then computed using the values derived: 
Scale Values = Grand Average - Row Average x Standard Deviation. 


These scale values are then transformed using a linear transformation. 


G. STEP BY STEP INTERVAL SCALE DEVELOPMENT 

The development of the interval scale as presented in [Ref. 7] is outlined below. 
For the present it is assumed that all p,,; are such that 0.02 <p,,< 0.98, so that the 
Z array is complete. The array P 1s defined as the cumulative frequencies of each. in- 
stance and z,, is defined as the values determined from a table of the normal distrib- 


ution for the p,, values. 


1. Arrange the raw frequency data in a table where the rows are scenarios (instances) 
and the columns are categories. Columns should be tn rank order with the least 
favorable category in the left colunin and the most favorable in the right colunnn. 


2. Compute the relative cumulative frequenctes for each row and record these values 
In a new table. This table is referred to as the P array and all values of p,,;> 0.98 
and p,,< 0.02 are removed. This creates an nx(m—k) array, where k is the 
number of columns removed. 


3. Treating these p,,; values as leftward areas under a Normal (0,1) curve, look up 
the values of Z from a table of the normal distribution. Record these as a new table 
Which will be the Z = || z,, || array for the computations that follow. 


4. For each instance, 1, in the Z array, compute the row average, 7, 


seeevor cach) columi j im the Z array, compute the coluinn average, b,. Note that 
b, is the value of the upper bound of category j on the scale being developed. 


6. Compute the grand average, b, of all values of the Z array. 


7 Compute B= > (b,— b)? , the sum of squared column differences. 


m—k 
j=1 


Hs 


m—k 
8. For each row compute 4, = ¥)(z,, — Zz)? , the sum of squared individual differences. 


J=) 
9. For cach scenane compute / = , an estimate of VG +c. 
10. Finally for each row (scenario) compute S, = b — Z, x |=. , for aller 


These S, ‘s are the scale values of the instances, on the same interval scale as the 
category bounds, b,. Now all instances and category bounds are on the desired scale 
and any linear transformation Y=a+fx NX, £>0, may now be performed to adjust 
the scale as desired, using the same transformation to move both scenarios and category 
bounds. [Ref. 7: pp.14-15.] 

In order to better understand this concept a short exainple will be presented. Raw 
data is first gathered and set up in matrix form showing the frequency that a category 
was chosen. The least desirable category appears in the leftmost column. For this ex- 


ample itis column A, and for the study conducted the letters correspond to the follow- 


ing: 
A = Totally Ineffective, 
B = Ineffective, 
C = Marginal, 
De sieeve 
E = Totally Effective. 


The first part of step 1 above is to record the raw frequencies (Table 1). The numbers 
in the left column correspond to actual numbers assigned to an instance of the pilot 


study. The pilot study is discussed in Section I below. 


Table 1. RAW FREQUENCIES 





The second part of step 1 1s to divide the frequencies by the total number of judges (20) 


who answered the questions to determine the relative frequencies (Table 2). 


mane 2. RELATIVE FREQUENCIES 





Step 2 involves deterinining the cumulative relative frequency of the array by summing 
each column with the values of the columns to the left. This table ts called the P array 


and ts given in Table 3. 


Table 3. CUMULATIVE RELATIVE FREQUENCIES 





os | 07 | 95 


The second part of step 2 is to remove all values of p,,> 0.98 and p,,< 0.02. Therefore, 
the column of 0’s and the column of I's is removed from Table 3, resulting in Table 
4, 


Table 4. REMOVE ZEROS AND ONES 





In step 3 the p,,s are treated as leftward areas under a Normal (0,1) curve. The values 
of Z, taken from the table of the standard normal distribution, are recorded as a new 
table of the Z array (Table 5). 


Table 5. ©.NORMALIZED VALUES 


Nc [8 | ¢ 
52 -1.64521 -0.841457 0.674189 


652 
= 652 






Step 4 involves computing the row average, z,, for each scenario, 1. 


In step 5, for each category j, compute the column average, b, , and note that 4, is the 


value of the upper bound of category j on the scale. 


Step 6 consists of computing the grand average 6. The results of these three steps are 


shown in Table 6. 


Table 6. NORMALIZED TABLE | 


No | Be eee ) 


0. 

eo 0.841457 ().674189 0.604159 
58 - 1.64521 -0.674189 1.64521 0.224729 
-1,64521 (0.674189 1.64521 0.224729 


1.64521 0.524002 0.174667 
-1.64521 0.524002 1.64521 0.174667 
eee ee) 


ee ee 


-1.64521 0.041309 [45101 bh = ~0.050965 





moh 


Instep 7, B=" — b)? , the sum of square column differences is computed as follows: 
j=l 


B =( —1.64521 —( —0.050965))? + (0.041309 — ( —0.050965))? + (1.4510! — ( —0.050965))" , 


B = 4.800061 . 


m— Kh 
In step 8 compute 4,= S\(z,,—2Z,)? , the sum of square individual differences (Table 7). 
j=l 


Table 7. SUM OF SQUARE DIFFERENCES 


1.08379 (1.056 3099 
ees) 2.01777 ~— | ~——«0.20201 3.49668 
| 66 | (3.49608 0.202014 
3 





In step 9 an estimate of the standard deviation \/a?+c is computed by ora 
“4; 


Finally in step 10, for each row (scenario) compute S, = b — Z x (= , for alli (Table 
2 | 
Table 8. SCALE VALUES OF INSTANCES 


NO. 


4. SO06056 


| a2 Y.744226 = —0.050965 — (—0.604159 x > 77497 ) 


4, SUELO 


58 my ee Es = gO ee 
= 0.155094 9.050965 — ( —0.224729 x 3 71646 





4 806056 


(Ns 
pee 5.71646 ? 


l 








—(.050965 — (0.224729 x goa 


a Foencos = _¢ 4, SO60S6 
0.212828 = 0.050965 — (0.174667 x ./=SCOPR® | 
92 -0.212828 = —0.050965 — (0.174667 x ae) 
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The column averages that correspond to the upper bound of the instances are given i 
Table > 


Table 9. COLUMN AVERAGES 
-1.64521 0.0413094 1.4510] 










A linear transformation Y=a+ Bx, B>0, may now be performed to adjust the 
scale as desired, using the same transformation to move both scenarios and category 
bends. In this example the transformation of the S,s was done by arbitrarily assigning 
a value of 100 as the upper bound of category D and zero (0) as the upper bound on 
category B. This results in the solving of two simultaneous equations with two un- 


knowns. These two simultaneous equations are: 


100 


ao + Bx (1.45101), 


() a + Bx (1.64521). 


Utilizing this transformation results in the values given in Table 10, 


Table 10, TRANSFORMED VALUES 
B C | D 
















46.262 
46.262 


The transformed values of both the category upper bounds and the scale values of in- 
stances are now all on one interval scale as seen in Fig. 2. This completes the transfor- 


mation of categorical responses to an interval scale. 


fle 





Marginal Effective 
0 54.47 100 


Figure 2. Transformed Values on Interval Scale 


Iti PILOT STUDY 

The judgemental! data base in this study involved the design of a (PILOT) questionnaire 
and its distribution to a representative group of Army officers at the Naval Postgraduate 
School. These Army officers will be referred to as “judges” in this study. The pilot study 
is used to determine the reasonableness and validity of the data gathering techniques. 
The study also allows changes to be made to correct any problems that may arise before 
sending the questionnaire out Army wide. The judges were asked to assume the role of 
the decision-maker of a Mechanized Infantry battalion with one mission: deliberate de- 
fense. It also provided the judges with the present STATE of the unit and the situation. 
The present STATE was reflected as a percentage of personnel, ammunition, combat 
vehicles and POL on hand or available. Two of these indicators (Personne! and Combat 
Vehicles) were evaluated at three levels (100%, 75%, 50%), and two (Anununition and 
POL) were evaluated at four levels (100%, 75%, 50%, 25%). This approach required 
3°) x (47) = 144 unique scenarios, each representing a different unit capability. 

In this pilot study 75 questionnaires were sent to U.S. Army officers (judges) with 
each questionnaire containing 48 different scenarios. A sample questionnaire is pre- 
sented in Appendix A. The judges were asked to determine, based on the situation and 
his units’ present STATE, to what degree his unit could accomplish its assigned mission. 
Irom each of the scenarios the judges determine, based on the STATE of his unit, the 
ability of the unit to accomplish its assigned mission. This ability was then rated in one 
of the following categories: 

ome Otay Ineilective, 


e Ineffective, 


e Marginally Effective, 
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e Effective, 


© foetaliy eileetive: 


In addition to the questions on cffeetiveness the judges were asked career questions, 
(e.g., branch of service, rank, staff experience), to determine the experience levels of the 
judges. The judges were also asked if they considered 50% for personnel and coinbat 
vehicles to be a lower bound for a unit being capable of accomplishing its mussion. 
These responses are presented in Appendix B, as a guide for future studies of this type. 
Of the 75 questionnaires 60 were returned and 15 were not returned until after the 
analysis had been completed. The analysis was conducted and the model constructed 
based on the first 60 questionnuires returned. 

[In general, the questionnaires required approximately 45 munutes to I hour to 
complete. Judges were asked to determine at what level their unit could accomplish its 
assigned combat mission based on the various combat situations presented. The exper- 
iment consisted of 144 different scenarios, and since the number of Army officers at the 
Naval Postgraduate School was limited, the experiment was organized into three groups 
of 48 scenarios cach. The scenarios that each group received and the order of the sce- 
narios in each group were randomly dctermimed with the use of a random nuniber gen- 
crator. This was done to avoid a trend in any given group, since it is undesirable for one 
croup to have a predominance of variables with high or low percentages. Each seenario 
had the same mission, eiiemy force and situation while the four independent variables 
were varied. In so doing, the relationslups of the independent variables and best fitting 


model was determuned as discussed in Chapter [V. 


I. TRANSFORMATION OF THE QUESTIONNAIRES 

Utilizing the procedure outlined in Section G, first the data was recorded and raw 
frequencies Were determined. The relative frequencies and cumulative relative frequen- 
cles Were then computed and are presented in Appendix C, page 54. With the eumula- 
tive rclative frequencies established the most difficult part of the procedure was to 
establish the groups of “like” categories (i.e., instances with responses only in categories 
A, B, and C are grouped with all other instances with responses only in categories A, B 
and C). All values of p,,2 0.98 and p,,= 0.02 were yemoved (1.e., all of the 1 sia 
O’s are removed in this study). These values are given in Appendix C. 

The next step was to normalize the values and transform them to one scale. This 
was accomplished using A Programming Language (APL) function that converts the 


cumulative relative frequencies to a scale value (steps 3 thru 10 as previously discussed). 


Ineffective Marginal Effective 
wo SL 50 WD. Shs 100 


Figure 3.  Jransformed Category Bounds on Interval Scale 


This APL function (called NORM) is given in Appendix I and the scale values of in- 
stances are shown in Appendix C. The final step was to transform all of the instances 
to the same scale. Another APL function (called TRANS) accomplished this transfor- 
mation of each individual group. To use this transformation two upper category bounds 
had to be established, the values of these bounds were arbitrarily given as 100 for the 
upper category bound for category D and 50 as the upper category bound for category 
B. The group with the categories B, C and D, was transformed first using the transfor- 
mation program producing a value of the upper categorv bound for C. These values 
were then used to solve for the values of the category upper bounds of the other groups. 
This puts all of the instances on one interval scale with the category upper bounds on 
the same scale. The transformed values are presented in Appendix C, page 63, and the 
corresponding category upper bounds are: 

e A (Totally Ineffective) = 33.32, 

e 8B (Ineffective) = 50.00 , 

* C(Marginal) = 72.933, 

e D (Effective) = 100.00. 


As seen in Fig. 3 the category upper bounds are now plotted on their interval scale. The 
transformed values of unit effectiveness were then used to conduct the analysis and de- 


termine the function that best fit the data. The analysis is discussed in Chapter IV. 


IV. THE ANALYSIS 


The preceding chapter explained the techniques used to determine the values of the 
dependent variable, unit effectiveness. This chapter develops the formulation of the 
model and shows how the model is used to support the Generalized Value System 


(GVS). In addition, an analysis of how well the model fits the actual data 1s explored. 


A. REGRESSION TECHNIQUE 
Once the transformed values were compiled, an analysis of the data and a deternu- 

nation of the best fitting function of the state variables and the values of unit effec- 
tiveness was done with the use of a regression technique in an A Programming 
Language (APL) package. Several functions that were tested using this regression 
technique: 

® Exponential, 

¢ Logarithmic, and 


@ Square root. 


None of these functions were determined to fit the values of unit effectiveness verv well 
since the standard deviation was unusually large. Appendix E shows the results of the 
functions tested and the coeflicient of determination (R Square) for each. An example 
of the logarithmic function is presented in Fig. 4. 

The next step in finding the best fitting function was to determine the shape of the 
four dimensional vector of personnel, ammunition, vehicles, and POL. This was done 
with the use of GRAFSTAT on the main frame computer at the Naval Postgraduate 
School. Each variable was compared to all other variables by plotting two variables 
against each other while holding the other two variables constant, using contour surface 
plots. These plots are presented in Appendix D. From these plots the general shape of 


the surfaces was determined to be elliptical; thus. the model form was given by: 
Yea +b (PER=A)Y toc CUO =O Faire Oe ee (4.1) 


Once the general shape of the surface was determined the next step was to perform 
regression analvsis on different forms of these functions to determine the best fitting 
model. The best fitting model as determined by the regression software (Appendix I) 


and utilizing the function in Eqn. 4.1 1s: 


"= 88.978 — .0056 x X1 — .0055 x 2 — .0054 x X13 — .0005 «x M4 , (4.2) 


where: 
c= 
X1 = (PER — 100, 

= (AMMO — 100), 

X3 = (VEH — 100), and 

X4 = (POL — 100). 


The values of unit effectiveness, 


~ 
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Figure 4. 
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Equation 4.2 has a coefficient of determination (R square) of 84% and a standard devi- 


ation of 6.41. [It was determined to be the best model among the techniques tested that 


could be developed from the data base. Dividing Eyn. 4.2 by 88.978 would set the unit 


effectiveness on a scale from 0 to l. For the purposes of this thesis this function (Eqn. 


4.3) is referred to as the Degraded Power Function (DPF), which ts the degradation of 


power based on the STATE of the entity at time, t: 
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aid 88.978 : oe) 


The DPF is now an estirnate of the overall unit effectiveness based on the four inde- 


pendent variables, personnel, ammunition, vehicles and POL. 


B. WILCOXON SIGN RANK TEST OF LOCATION FOR PAIRED SAMPLES 
Based on Eqn. 4.2 a Wilcoxon Signed Rank test was used to determine how well the 

model fit the actual transformed values. For this test X represents the transformed val- 
ues and Y represents the predicted values using Eqn. 4.2. Since the direction of the 
difference between X and Y could not be anticipated, the hypothesis is set with a two 
sided alternative for the median of D = X - Y. and it is assumed that the differences, 
D, are symmetrically distributed about the median of D ( Jf, ). Therefore, the null- 
hypothesis (H), and alternative hypothesis (A), are: 

ht = Ulead 

A: Vi, = 0: 


Using the Wilcoxon Sign Rank test, a table was established to determine thie dif- 
ferences between the Transformed values (X), and the predicted data (Y). The sign of 
the difference between X and Y is used to determine the sign (positive or negative) of the 
rank of the D's = X,— Y, for all instances1. The ten steps in this process are explained 


below. 
1. Find the differences between the actual data and the predicted data. 
2. Find the absolute value of these differences, noted as | D | . 


3. Rank the | D,| with the smallest absolute difference having a value of I, the next 
smallest a value of 2, etc. 


4. Give the sign of the D, to these ranks. 

5. Find the value of T, = >& positive ranks. 

6. Find the value of 7_ = § negative ranks. 

7. Forn = 15 determine the larger valve ol eam 


Ss. For n> 15 determine the larger value of Z,, and Z_, using Eqn. 4.3 below jean 
use this value for a right-tail probability look up for the p-value, sometimes called 
the eritiedisiever 


9. Using Table G or A (when using normal approximation), (Ref. 8: pp.123-[4]] 
determine the p-value in terms of a nght-tail probability for one of ( 7, , T_ ), or 
CANA alec! 1 Sly) 
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Minis f-valwe 15 the basis Of acceptance or rejection of the Null Hvpothesis (1), 
which stated the median differences between the transforined values (X) and pre- 
eireted { 1 Jeralues are Zero. 

Based on the Wilcoxon Signed Rank test, presented in Appendix G with n = [44 the 
Pimwot tanks were: J. = 5397 and 7, = 5043, aud Z,, and Z_, were computed using 
Eqn. 4.3 


Fis a 7 -05- SU 
Me lg (43) 
+,R na + Don + 1) ue nia + W202 + 1) 
24 24 


which resulted in Z.p = -0.354 and Z_, = 0.352, Z_, was used since it had the larger 
value. Based on the information !n a normal CDF table, and since this was a two-sided 
hypothesis test, a p-value of 2x(0.3624) = 0.7248 was computed. Based on this p- 
value, the null hypothesis could not be rejected and the assumption that the predicted 
values approximate the actual values was valid. Note that the most deviation was at the 
large values of X (Appendix G). Since this model will not be used to predict values of 


the state variables above 100%, this model is considered quite accurate. 


C. DEVELOPMENT OF A MAPPING FUNCTION 
Recall from Chapter IJ that Predicted Adjusted Basic Inherent Power (PABIP) may 
be calculated using Eqn. 4.4 


PABIP(X\( 1 )|SX1( 1, )) = ABIP(SX\(t,)) cme Gy pafor Tt = ie. (4.4) 


In order to demonstrate how Eqn. 4.4 relates to time, a short example is used to de- 
scribe how PABIP functions. If the unit has the following constants, the shape that 


Equation 4.4 takes on 1s depicted in Fig. 5, 


e L = 0.02, a rate of change which may have several components such as attrition, 
reinforcements, and logistics effects, 


fee, = OU, the end of the planning horizon, 


° ¢, = 0, the point in time that a unit is in position to accomplish the mission, 


cae | Ue tne Current time, and 


e ABIP = 1000, 


Note that PABIP( t) = ABIP, when t=0 and the formula for PABIP (Eqn. 4.4) is 


only based on the value of ABIP at one specific time (t = ¢,). For 1>4, Eqn. 4.4 uses 
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an exponential decay to forcast PABIP. Therefore, the state of the unit 1s only used at 


ume, ¢,. 


The proposed model to replace Eqn. 4.4 1s: 
PABIP(XM(t)|SX1(t)) = ABIP(SX1(t)) x DPE(SX1(1)), (4.5) 


thereby eliminating the exponential decay function, e-'-%) , and utilizing the state of 


the entity at discrete incremenisiah wimle, t. 1On mi aie 


PABIP USING EXPONENTIAL DECAY 
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Figure 5. —‘ Predicted Basic Inherent Power using Eqn. 4.4 


Also from Chapter II the Situational Inherent Power (SIP) is computed by using 


Equation 4.6 
SIP(XI(£)\(t)) = PABIP(SXI(1)|SK1(t,)) xe for OS t Sty, (4.6) 


SIP(X1(t)) = PABIP(SXI(t)) for t>t,. 
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Recall that the rate at which a unit attains readiness (mission position discount factor), 
D, accounts for the tiine until the unit is in position to perform its mtssion. The pro- 
posed model could make the SIP more accurate since it continually uses the forecasted 
fiaivero! the unit to calculate the PABIP. Essentially the PABIP will equal ABIP if the 
forecasted state of the unit is continually updated, as ALARM has the capability of 
doing. Thus can be seen in the following example. 

A Blue forces is defending against an attacking Red force. The Blue force has a 
Basic Inherent Power (BIP) of 1000 and the Red force has a BIP of 2000. Both forces 
Pewee 100% streneth with the Red force 16,200 meters from the Blue force at t=0. 
fate o! inarch for Red force is 270 ineters per minute, therefore, the Red force is pres- 
ently 60 minutes from the Blue force’s position. In order to use the SIP, D (mission 
position discount factor) must be coinputed based on the time it takes the Red forces to 
reach the position to conduct an attack against the Blue force (for this example this tune 
will be 30 minutes). A negligable amount of power is assigned to the Red force at t=0 
(i.e., 5% of the Red force PABIP). This grows until the Red force is at 100% of PABIP, 
when they are 30 minutes from the Blue force’s position. Thus the discount factor for 
the Blue force considering the Red force is determined by solving 0.05 = e-%) for 
D, resulting in D = 0.0998 / minute. Since the Blue force 1s in position, their STP(t) 
= PABIP(t) for all t. During the conduct of this battle tlle Blue and Red force’s SIP’s 
will be plotted, as well as the difference between their respective SIP’s. For the first 30 
minutes of the battle the attrition was due mainly to artillery and vehicle break-downs. 
The rates of attrition for Blue and Red forces are presented in Table 11. Based on the 
famesesivown in lable 11, the values for Blue and Red force's DPF ’s, PABIP’s, SIP’s. and 
the difference between the Blue and Red forces SIP’s were computed. These values are 
presented in Appendix H, and a graphical representation of the battle is shown in Fig. 
6. 


Table Il. BLUE AND RED FORCES ATTRITION RATES PER MINUTE 


Red Forces 
( Viin) PER | AMMO] VEII PER | AMMO! VEII POL 
s-.30 0.000 0.167 0,067 











Time 
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Note that as the Red force moves into attack position their SIP increases but their SIP 


is also being degraded by the attrition factors. 


BLUE VERSUS RED 
SIP AS A FUNCTION OF TIME 
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rigure 6. Blue versus Red Battle 


Once the Red force begins the attack, the attrition rate increases and their SIP re- 
duces at a greater rate than the Blue force. However, at t= 60 minutes the Red force's 
SIP is greater than the Blue force’s SIP, which would indicate that the Red force would 
win the battle. 

If the Blue commander could forecast the state of both Red and Blue forces at 
t= 60, he could request that reinforcements arrive prior to the Red attack. The same 
eXxainple is run again with the same attrition rates, presented in Table 11, and the same 
BIP’s for both forces, only this time the Blue force requests that reinforcements arrive 
by t=30 minutes. The reinforcements in this example have a BIP of 1000, Tlus second 
battle is presented in Fig. 7. Vhe Red force’s SIP is increases until they are in position 


to attack and the Blue forces SIP is increasing as the reinforcements move into position 
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BLUE VERSUS RED 
SIP AS A FUNCTION OF TIME 
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Figure 7. Blue Versus Red Battle with Reinforcements 


to support the Blue foree. At t=30 nunutes, both force’s SIP’s are almost the same. 
As the Red force attacks, their attrition rate 1s higher than the Blue force, enabling the 
Blue force to win this battle. Therefore, if a decision-maker can furecast the anticipated 
State vector values of the Red and Blue forces over 4,5 ¢<¢,, he can determine by what 
time resupply or support has to occur in order to positively affect the outcome of a 
battle. 

This example illustrates some very important aspects of representing and modeling 
the decision process. First, future state forecasting is the only way the Blue force could 
decide to initiate commitment of the reserve at t=0 to realize success at t=60. That is, 
decision tables based on state vector values at t=V could not cause this decision to oc- 
cur. Secondly, the forecast is based on anticipated force states for 4. <¢<¢; using the 


DPF and not on an arbitrary power rate of change. Finally, this process is very close 


to the process actuallv used by the $-2 (INTEL), S-3 (OPNS) stafF functions in an actual 
combat unit. These techniques may be usable not only as an aid in systenuc combat 


simulations but as a decision aid to a unit staff. 


D. FUTURE APPLICATIONS 

The feasibilitv of developing a mapping function between a multiattribute state 
vector and a unit's power degradation has been demonstrated using a Pilot experiment. 
This methodology has potential use both in combat models and as a decision aid to 
combat unit staff planners. 

It is recommended that this same experiment be repeated using subjects currently 
in mechanized battalions to determine the validity between that group and the Pilot 
group. Next, additional experiments must be designed and implemented for larger units, 
different missions, different threat forces, and possibly a revised set of ind pendent var- 
lables. For example, considerations of a division will likely require additional inde- 
pendent variables such as Class IX (repair parts). Environmental factors such as 
weather may also be included. 

Obviously the goal is to accomplish the required mapping with a minimum number 
of independent variables. Ilopefully, this thesis provides the basic methodology for 
conduct and analysis of a large famuly of experiments required to realize robust mapping 


functions. 


30 


APPENDIX A. SAMPLE QUESTIONAIRE 


A. PURPOSE & MOTIVATION 
The purpose of this questionaire is to obtain an estimate of the degradation in a 
units effectiveness based on the threat posed, by the degradation in the 4 kev variables 
of a unit. This questionaire will not just be used as a subject for a thesis and later 
diregarded. It is to be used to help the S-3 and decision maker to deternune if he can 
handle a mission based on 4 key factors. Your answers are therefore very important to 
insure a good decision is made. You answers, as the decision maker, of the percent de- 
gradation of your units effectiveness based on the changes to these Key variables 1s a 
measure of the relative importance of each variable to the accomphshment of your 
mission. Your answers will also help to develop a more accurate and realistic represen- 
tation of how changes to the key variables effect your view of its relative importance. 
I. Key Variables 
The 4 kev variables used throughout this questionaire are: 
“> Personnel, 
O% Ammo, 
% Weapon Systems / Combat Vehicles, and 


poe OL (Fuel) . 


2. Instructions 

In the remainder of this questionaire, you will be asked to place yourself in the 
role of the decision maker of the unit and determine how changes to your units fucl, 
ainmo, personnel and vehicles will effect your interpretation of your units ability to ac- 
complish its assigned mission. 

Please respond to the questions asked tn accordance with your feelings regard- 
ing the situation. There is no nght or wrong answer to any of the questions. As a de- 
cision maker in a combat situation you will be required to make rapid estimates of the 
situation. Therefore, with this in mind vou should only take enough time to fully un- 
derstand the situation presented and record your response. Once you have recorded a 


response you should not change your response. 
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You will receive an example question and two practice questions. You are then 
asked to answer 48 questions by putting an X under the response you feel best de- 


scribes the unit’s ability to accomplish its assigned mussion. 


Based on the following situation you will be asked to answer questions on the 
degree to which you determine vour unit is able to accomplish its assigned mission. Use 


the following situation to answer all of the questions in this survey. 


Situation 


1. Enemy - 127th Motorized Rifle Regiment. 


2. Friendly - Your unit 2nd Bn 41st Inf Mech is currently conducting deliberate de- 
fensive operations along the forward line of troops (FLOT). 


e Your unit is presently in prepared defensive positions. 
e Your unit is the forward unit. i.e. no units to vour front. 


e Your unit is currently engaged in combat with the enemy. 


Nlission 


Your unit 2/41st Inf (M) will conduct a deliberate defense of present positions 
for a nunimum of 24 hours, longer if possible, to prevent the enemy from controlling this 


key terrain. 


Based on the above scenario and mission answer the following questions. 


a2 


B. EXAMPLE OF THE RESPONSE FORM 


socdedescatcrcakcdeseveseskeslese seve dcsdedestvdevevesevesvedevedescscdescdedesnevesdedescdevedcdentdededescdedoucdesescdededcdenesedescdevedouesedve 


THE CURRENT STATUS OF YOUR UNIT IS : 


75% PERSONNEL, 
25% AMMUNITION, 
50% WEAPON SYSTEMS, and 
50% POL (FUEL) . 
BASED ON THIS STATUS, INDICATE BELOW THE CURRENT EFFECTIVENESS 


OF YOUR UNIT'S ABILITY TO CONTINUE TO ACCOMPLISH ITS CURRENT 
MISSION OF DELIBERATE DEFENSE. 


TOTALLY TOTALLY 
EFFECTIVE EFFECTIVE MARGINAL INEFFECTIVE INEFFECTIVE 


devededevesedededcdenvcdevedevedededcdedekdeddcdevesevedcteocdcdcdedcdedcdedcdededcdedededvcacdedcddesdcscdevesovlevevescdedcdedede dededcvede 
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APPENDIX B. PERSONAL HISTORY FORM 


Please complete the following form by entering the appropriate responses: 

1) Time on active duty (average): (8) Years, (9) Months. 

2) Time spent as a staff officer (average): (3) Years, (2) Months. 

3) Present Rank ( O=LI@ (7A eG so er le: 

4) Were the situations presented understandable? ( 55 - Yes), (5 - No). 

5) Would you say the situations were realistic? ( 49 - Yes), ( [Ll - No). 

6) You may have noticed that the Personnel and Weapon systems are only evalu- 
ated at 100%, 75%, and 50%. Would you consider the 50% a minimum level to eval- 
uate these two variables or should one or both be evaluated at 25% ? 

(check the statement that applies). 

(2) Just Personne! should be evaluated at 25%. 
(2) Just Weapon systems should be evaluated at 25%. 


(10) Both Personnel and Weapon systems should be evaluated at 25%. 


(46) I consider 50% for both a minimum level. 
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APPENDIX C. TRANSFORMING CATEGORICAL RESPONSES TO 


INTERVAL SCALE 


Recall from Chapter III the procedure for Constructing Interval Scales from 


Categorical Responses 1s: 


l. 


2B 


Arrange the raw frequency data in a table where the rows are scenarios (instances) 
and the colunins are categories. Columns should be in rank order with the least 
favorable category in the left column and the most favorable in the right column. 


Compute the relative cumulative frequencies for each row and record these values 
In a new table. This table is reffered to as the P array and all values of p,; = 0.98 
and p,;< 0.02 are removed. This creates an m xX (#1—k) array, where k is the 
number of columns removed. 


Treating these p,,; values as leftward areas under a Normal (0,1) curve, look up 
the values of Z from a table of the normal distribution. Record these as a new table 
which will be the Z = || z,, || array for the computations that follow. 


or each instance, 1, im the Z array, compute the row average, z, 


For each column j in the Z array, compute the column average, 6,. Note that 
6, 1s the value of the upper bound of category j on the scale being developed. 


Comipute the grand average, b, of all values of the Z array. 
mn—k 


Somputce B= > (b= b)? , the sum of squared column differences. 
Me m—k 
For each row compute 4,= >) 3) (z,;—z,)? , the sum of squared individual differ- 
ences. tes 


B ieee 
For each scenario compute 7 > an estimate of ja?+c¢ 
“1, 


10. Finally for each row (scenario) compute S, = b -—Z,x zs , for all 1. 


A 


{ 


iiiccems. S are tiie scale valuies Of the Instances, On the same interval scale as the 


category bounds, 5,. Now all instances and category bounds are on the desired scale 


and any linear transformation Y=a+fx X, fB>0O, may now be performed to adjust 


the scale as desired, using the same transformation to move both scenarios and category 
bounds. [Ref. 7: pp. 14-15.] 


The remainder of this appendix shows the steps used to transform the categorical 


responses to an interval scale. 


NS. 


GROUP | RAW FREQUENCIES 


A. 
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APPENDIX D. PAIRWISE INDEPENDENT VARIABLE COMPARISONS 


The purpose of this appendix ts to show that plotting each pair of independent 
variables results in elliptucal shaped contour lines when the other two independent vari- 
ables are held constant. The case of Ammunition plotted with Vehicles 1s explained be- 
low. The x and v axes refer to the independent variables being plotted while holding the 
other two independent variables constant. In Table 12 Ammunition 1s compared with 
Vehicles while holding Personnel and POL at 75%>. The dependent variable’s value 1s 
shown as a function of the independent variables in Fig. 8. 

mom this 3-ID plot Fig. 9 reflects the contour plot where the cuts are made at de- 
Pemdent variable values of 50.60,70,80, and 90. These cuts are labeled in Fig. 9. The 


reinaining pairs of independent variables are shown tn Figures 10-21. 


Table 12. DEPENDENT VARIABLE AS A FUNCTION OF ANIMO AND VET 
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Figure 8.  3-D Contour Plot 
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Figure 9. Ammunition Versus Vehicles Holding Others at 73% 
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Figure 10. Personnel Versus Ammunition Holding Others at 100% 
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Figure 11. Personnel Versus Ammunition Holding Others at 75% 
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Figure 13. 
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Figure 14. Personnel Versus POL Holding Others at 100°% 
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Figure 16. Ammunition Versus Vehicles Holding Others at 100% 
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Figure 17. Ammunition Versus Vehicles Holding Others at 73% 
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Figure 18. Ammunition Versus POL Holding Others at 100° 
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Figure 19.) Ammunition Versus POL Holding Others at 73°o 
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Figure 20. Vehicles Versus POL Holding Others at 100% 
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Figure 21. Vehicles Versus POL Holding Others at 75% 


APPENDIX EE. REGRESSION TECHNIQUES 


The results of regression of the values of unit effectiveness (YY) and the independ- 
ent Variables (PER, AMMO, VEH, POL) are shown tn this appendix m an attempt to 
determune the best fitting function. 
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APPENDIX F. 


Regression of the values of unit effectiveness (YY) and the 


BEST FIT REGRESSION MODEL 


independent variables 


(PER, AMMO, VEH, POL) was conducted resulting in the following best fitting model. 


X1 = (PER — 100)?, X2 = (4MMO— 100) 
X3 = (VEH — 100), X4 = (POL — 100) 
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APPENDIX G. WiILCOXON SIGN RANK TEST 


The transformed values are represented by X and the predicted values are repres- 
fmiedspy Y. The Wilcoxon Sign Rank Test for Location of Paired Samples tests the 


hypothesis that the median difference between X and Y ts zero. 
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APPENDIX H. BLUE VERUS RED BATTLE 


This appendix presents a battle between Blue and Red forces. Section A is a battle 
with no resupply or reinforcements, and Section B is a battle with the Blue receiving re- 


inforcements by time, t= 30 minutes. 


A. BLUE VERSUS RED BATTLE WITH NO RESUPPLY 


Time Blue Red Blue Red Blue Red (Red - Blue) 
(min) DPF DPF PABIP PABIP Sue SIP SEP s 
0 1. OLIG 1.0000 1000 1500 1000 100 -900.0 
2 1. 0000 IES SOLO LO O07, 1500 1000 122. 1 Vase 
4 0. 9999 L200 GG 959529 1500 999.9 149.1 -850.8 
6 0.9997 0. 9999 999.7 1500 9997 LS2.0 Sole 
8 0.9995 C7998 999.5 1500 STS S 222.3 ee? 
10 0.9992 0.9997 999.2 1500 999.2 271i. 4 Pea 
12 0.9988 0.9996 998.8 1499 998.8 331.3 "=s6675 
14 O2o9es 0.9994 998.3 1499 998.4 404.5 =593.9 
16 059979 029992 997.8 1499 997.9 493.8 <504.1 
18 0.9973 0. 9990 997. 3. 1499 997.3 602.8 -394.4 
20 0.9967 OL SI Tals 996.6 1498 996.7 735.9 26057 
ps 0.9949 0.9956 994°-9 1498 994.9 893.4 -96.5 
24 0.99278 0.9983 DO ao 7 992.3 WHeeo? . 0 Oe 
26 0.9888 0.9980 988.9 1497 988.8 1339.0 349.9 
28 0. 9844 Og] DS ele oG 984.4 1634.0 649.7 
30 0.9791 029974 Se Sl By IEE 979. ie 1955 10 62-0 
Be 0.9573 0.9901 968.9 1482 96773 oe POO 
34 ORR) Ss at 0.9768 Gao), 1457 952.1) es4 1TOO2. 0 
eG 0.9335 0.9576 G30 o, 1420 933. 509 n5 981.7 
38 Q.. 95 0.9324 G20 be 1372 911, Selo Soc ane 
40 0. 8862 OF “20s S98.7) 132 886.2 1803 ONG: 
i O57 5 0. 8642 874 1240 857.5 Zs 67 029 
44 One2 55 One212 BiG, Sool lon 825.5 mero 816.9 
46 0.7901 Q.77e22 816 1062 790.1 1544 FS as: 
48 Ones Ge vel7 3 OZ 955 751.3 Sela BSSe2 
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ae O66 37 0.5895 700 7060.7) 9660547 eee? oy lo ate 
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56 0.5626 0.4379 62028 Pes) Ge, 8) Seas Sore 2 
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60 0.448 O° 2624 52256 70.4) 4a eae Hie ars: 


B. BLUE VERSUS RED BATTLE WITH REINFORCEMENTS 
Blue forces receive reinforcements by time, t= 30 minutes which is reflected by an 


imecrease in the Blue forces PABIP, and SIP. 


Time Blue Red Blue Red Blue Red (Red - Blue) 
(min) DPF DPF PABIP PABIP SIP SIP SIP's 
0 0000 sl. 0000. 1050 2000 #1050 100.0 -950.0 
2 1.0000 1.0000 1061 2000 £1061 122.1 -938.9 
4 0.9999 1.0000 1074 2000 £1074 149.1 -925.3 
6 0.9997. 0.9999 1091 2000 1091 182.0 -908.7 
8 0.9995 0.9998 1111 2000 tt] 222.3 -888. 3 
10 0.9992 0.9997 1135 1999 =1135 Dee. Carn 
12 0.9988 0.9996 1164 1999 1164 331.3 -833.0 
14 0.9984 0.9994 1200 1999 1200 204. 5s 7 5, 9 
16 0.9979 0.9992 1244 1998 1244 COGS S750 5 
18 0.9973 0.999 1298 1998 1298 602.8 -695.2 
20 0.9967 0.9988 1364 1998 1364 735.9 -627.7 
pO 0.9949 0.9986 1443 1997 1443 898.4 -544,. 3 
24 0.9923 0.9983 1539 1997 1539 1097.0 -441.9 
26 0.9888 0.998 1655 1996 =: 11655 1339.0 -316.4 
28 M0846 0.9977 1797 1995 iwoy 1634.0 -163.0 
30 0.9791 0.9974 1970 1995 1970 1995.0 -24.3 
32 0.9673 0.9901 1935 1980 81935 1980.0 Coa 
34 0.9521 0.9768 1904 1954 1904 1954.0 V0), 
36 0.9335 0.9576 1867 1915 1867 1915.0 48.0 
38 0.9115 0.9324 1823 1865 1823 1865.0 42.0 
40 0.8862 0.9013 1772 A10)S) ae yy 1803.0 31.0 
Dg) Oee575 0.8642 1715 ies, 1715 l723e0 ey 
44, OeePoomeeOe S212 1651 1642 1651 1642.0 -8.6 
46 Dey oO0b 80.7722 1580 1544 1580 1544, 0 =35.8 
48 Om75i3)) 0.7173 ° 1503 1435 1503 1435.0 -68.1 
50 0.7092 0.6563 1418 Shs aaales 1313.0 -105.7 
52 0.6637 0.5895 1327 ake ase: 1179.0 -148.4 
54 0.6148 0.5166 1230 LOS ele so 1033.0 -196.3 
56 Ooo 0.4379 1125 BS Pe BS 875.7 =-249.4 
58 Ooo 0.5531 101g yA). yuc wen ROR EA 706.3 -307.7 
60 0.4480 0.2624 896 524.9 896 S7Geg S37 lal 
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APPENDIX I. PROGRAMS DEVELOPED FOR THESIS 


This appendix contains the A Programming Language (APL) programs used in this 


thesis. The following programs are available at the Naval Postgraduate School. 


A. STANDARD NORMAL TABLE 


782 1.781477937  1.821255978 i va3s02/ 42 
)*0.5 


B. STANDARD NORMAL QUANTILE 


VNQUANCOIV 
V come P:A3:B8:C:D 
} +(€(+/(A<(2s0)v (P27 )) 20) /nd 
i Ce 2.515517 0.802853 0.010328 
1] De 1.432788 0.189269 0.001308 
] P<((A«(PS0.5))xP)+((P>0.5)x(1-P)) 
J B<+(@Px 2)*0.5 
oe (2x42 )-1)x-B-((Bo.* O 1 2)4.xC)#(1+°(0(B°.* 1 2 3)+.xD)) 
] E IS NO QUANTILE FOR P = ',8A/P 


on || | [| Ue | Ce (| | CC eC CU | |e | | || (|e (| oe |e (| oe | || ee (| en |e | | 
MOF EFEFE EE FE WWWWWWWWWWHONNDYNNDNDANNEP RRP BREE EO ONOU Ewe 
RPODOONOMNFEWNHROWOONOOEWNROWONOOMEWNROW ONO FWP OW Wott tty 


Lt itt 
+ 


The following programs were written by the author for the data transformation and 


data analvsis portions of this thesis. 


C. NORMALIZATION OF CUMULATIVE EREQUENGIES 


VNORMCOIV 

v WORN MTX;NOR:ROWAV:GRAV:S >AAI;II 
O+'7NPUT THE CUMULATIVE bbeaue hedest 
Het TN SURE THAT THERE ARE NQ VALUES OF' 
N<+'ZERQ (0) OR OF ONE (1)! 


MT X<«Q 
Eee MTX 


PORAVe (+/NOR)+($02]) 

ROWAV<(S(11,1)pROWAV 

BomAV~(+7NOF)+(SL1]) 
GRAV<(+/(+/NOR ))+( (S013 )x (S02) )) 

O<' NORMALIZED VALUES ROW AVERAGE? 


i+ 

DHF, »ROWAV 

f COLUMN AVERAGES ! 
< 


fie! SCALE VALUES = GRAND AVERAGE - (ROW AVERAGEX(BtATI)*. 3)! 


OSSST, GRAV ,ROWAV ,SQR 
pe! COLUMN AVERAGES | 
= 


D. TRANSFORMATION OF DATA 


VIRANSCOIV 
THANS;} UP ; LON SBETA | eee ies 5 OO a arc ma tad 


1] cous THE COLUEN AVERAGES FROM THE NORM FUNCTION'! 

33 Paes 

cu) UPeaGnnsa 

5] “OLLOWSCOLT A ‘le 

£8: Get in PUT THE VECTOR YOU WANT TRANSFORMED! 

8 Us | INPUT THE UPPER NUMEER YOU WANT WITH THIS GROUP! 

6194 Us! INPUT THE LOWER NUMSER YOU WANT WITH THIS GROUP! 
< 

ay BEPA<(UR=LGOW )+ (GOLUPSCOuLOu. 

[13] ALPHA<UP-BETAxCOLUP 

[14] TReALPHA+BETLAXHX 

[151 PRCOL<ALPHA+BETAXxCOL 

F416] G+ | TRANSFORMED COLUMN UPPER BOUNDS ARE NOW! 

f17] O+TRCOL 

r137 

peg 

[20] 


Ae: TRANSFORMED Ae 
O<7R 


90) 


OWONIAMNEWNROWOND UF Wh OMI 


WMNONNNNNNNNNRPRPRP PRP RPRRPREWOONOMFEWhHOrF 
ar | |] | | | WY | US| | OY WW] WY 


soe | gee | Gee | ee | | | fee | coe | cee | mee | cee | nee pee | coe! | co | gene | see | gee | gene | cee | moe | pee | oe | me | ee | cee | ee | ee | eee | 


WILCOXON SIGN RANK TEST 


WWILCOXCO]y 
V WILCOX: X%:Y:DIF;ABSDIF;RANK:SIRANK;SHP;SI;:SJ;SDY3;N321;22 
ake INPU? bnb’ x VALUES | 


q-! INPUT THE Y VALUES '! 
<< 

DIF+X-Y 

ABSDIF+|DIF 

SHP=0X 


RANK ha ABSDIF) 
RANK<(SHP )oRANK 
STRANR=RAN x (DIF +ABSDIF ) 


M<' Xx y D 

Ax Y,DIF,ABSDIF,RANK,SIRAN 
TMIN«+/SI<SIRANKx (SIRANKK<O ) 
TMIN+«+/|7M 
TPLUS<+/SJ<SIRANKx (SIRANK>O) 
TPLUS<+ TPLUS 

N<SHP 


( 
LUT N< ee )= (0. 


DIF | | RANK | [SIGN RANK|! 


1235528 +SDY 

£2~CAUSS ZPLUS 

et PhUs a shes ss; [T MINUS = ',(8IMIN) 
<2 PLUS =, (SZev0S) Z MINUS = ',(8ZMIN) 
aoa P- VALUE FROM Z PLUS = ',(8Z2) 

Oe: p- VALUE FROM Z MINUS = ',(82Z1) 


a! 


F. CALCULATION OF CUMULATIVE FREQUENCIES 


VIUDGESCOIY 
JUDGES;D:N;SHAPE: FREQ; 3s CUMEREQ 

het HOOT hire RAN DA OINT 

fe' HS MURB EAs Beery aPTOnS | 

O< Ti EACH ROW MUST BE THE SANE! 
< 


N<+/D 
ite see 


NeNC1 
FREO*«DLD:1(SHAPE[2])] 3 
Q< | RELATIVE FREQUENCY 


CUMFREQ<+\FREQ 
O<« ' CUMULATIVE FREQUENCY! 


(]e_l ono Seee eee == eee 
= CUNFREQ 
V 


[ae | oem | coe | cee | eee | cee | ee | cee | ee | ee | ee | ee | ee | ee | ee | ee | ee | 
F32R EE PPEEIQCONIOUN FWP 
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